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Predictions are made for the forward-backward and charge asymmetries in bottom-quark pair 
production at hadron colliders. Tree-level exchanges of electroweak (EW) gauge bosons dominate 
the Standard Model (SM) contribution to the asymmetry near the Z-pole. The mixed EW-QCD 
corrections are computed in an approximate way, and are found to be small in magnitude. These 
SM predictions are consistent with experimental results from CDF, DO, and LHCb. In particular, 
CDF and LHCb find that the asymmetry in the invariant mass bin containing the Z-pole is larger 
than in the adjacent bins, as predicted. Several beyond the Standard Model scenarios proposed for 
the top-quark forward-backward asymmetry, including a 100 CeV axigluon, are disfavored by this 
combination of SM predictions and measurements. On the other hand, modified Zhb couplings can 
explain the 2 g discrepancy in the bottom-quark forward-backward asymmetry at LEPl, while being 
consistent with the results of CDF and LHCb. It is also shown that f-channel W exchange makes 
a non-negligible contribution to the charm-quark charge asymmetry. 


I. INTRODUCTION 

For several years, measurements of, and related to, the 
forward-backward asymmetry in top-quark pair produc¬ 
tion (Apg) at the Tevatron were consistently higher than 
the Standard Model (SM) predictions [T}1S]. In particu¬ 
lar, the measurement of the forward-backward asymme¬ 
try at high invariant mass {Mtf > 450 GeV) was quoted 
to be 3.4 standard deviations above the next-to-leading 
order (NLO) QCD prediction [T]. This result sparked 
much work - both theoretical and experimental, within 
and beyond the SM (BSM) - trying to find the source of 
the discrepancy. See [5] for a review. 

One proposal to help solve the problem was to mea¬ 
sure the forward-backward asymmetry in bottom-quark 
pair production [311], and the analogous charge 

asymmetry at the LHC (Aq) [S], see also [10]. Such a 
measurement would likely be difficult due to dominance 
of gluon fusion initiated bottom pair production, which 
does not generate an asymmetry, over qq initiated pro¬ 
duction, the dominant top pair production mechanism. 
In addition, there were expected to be further complica¬ 
tions due to 6-tagging inefficiencies, and processes such as 
neutral H-meson oscillations and cascade decays spoiling 
the correlation between the charge of the decay prod¬ 
ucts and the charge of the parent 6-quark. However, the 
upside would be valuable information about the flavor 
structure of the source of the top asymmetry. Calcula¬ 
tions of the bottom-quark Afb in the SM including some 
electroweak (EW) effects had been made both before [TT] . 
and after [12] these proposals. 

In [13], it was realized that a Z-boson decaying to 
66 would have significant consequences for the analysis 
of the bottom-quark asymmetry. [63] It was found that 
near the Z-pole, tree level exchanges of electroweak gauge 
bosons dominated the SM contribution to the Afb- This 
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is unlike the top asymmetry where NLO QCD is the lead¬ 
ing contribution. Far enough above the Z-pole, NLO 
QCD does provide the leading SM contribution to the 
bottom asymmetry. 

The effects of various BSM scenarios on A’^p, and their 
relation to Apg, were also investigated in [l3].[64] How¬ 
ever, the tt asymmetry discrepancy has recently been re¬ 
solved through a combination of theoretical and experi¬ 
mental work, namely results from DO using the full Teva¬ 
tron Run-2 dataset [I4HT6] , and a SM prediction of A^'pg 
at next-to-next-to-leading order (NNLO) in QCD [17j . 
Nevertheless, in the interim, measurements of A^Jlg by 
the CDF |18l 119] and DO collaborations [101121] and of 
by the LHCb [HI HO] collaboration were made, and 
a new ~ 3cr discrepancy has appeared. 

The preliminary results of CDF at high invariant 
mass [18] are both consistent with zero, and agree with 
the SM predictions of [TO]. In addition, CDF was able 
to exclude a wide axigluon with a mass of 200 GeV as 
an explanation of A*pg, while not excluding an axigluon 
with a mass of 345 GeV. The preliminary results of CDF 
at low mass [TO] agree with the SM predictions of m , in¬ 
cluding a measured asymmetry in the bin containing the 
Z-pole that is larger than the asymmetry in the adjacent 
invariant mass bins. 

LHCb updated their preliminary analysis of 7 TeV 
data [22] after Ref. [13] came out to include an additional 
invariant mass bin centered on the .2-pole. In their pub¬ 
lished result. Ref. [20], they find that central value of 
in this 2-pole bin is the largest of the three bins in 
the analysis, similar to the prediction for A’^pp in m- 
However, there was no dedicated SM prediction for this 
measurement at the time it was released. 

DO did an analysis (201 El] of the forward-backward 
asymmetry in charged H-meson production, Afb {B^), 
rather than a jet based analysis like CDF and LHCb. 
Clearly, there will be none of the charge tagging system¬ 
atic issues discussed earlier in an analysis with charged 
mesons. Also, because DO used an inclusive sample of 
they have more events and thus a smaller statistical 
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uncertainty as well. However, this inclusiveness comes 
at the price of including far more gluon fusion initiated 
events, which dilute the already small asymmetry. DO’s 
result is (H±) = [-0.24±0.41(stat.)±0.19(syst.)]%. 
This is consistent with zero. However, this is also 
3.3 standard deviations below their SM prediction of 
Ap^ {B^) = [2.55 ± 0.76]%, which was made using 
MC@NLO+Herwig [53]|2S]. DO also finds that the mea¬ 
sured asymmetry is lower than their SM prediction for 
all pseudorapidities, and for px (B^) =9 — 30 GeV. 

DO has also made a measurement of the forward- 
backward asymmetry in A° and production |26j . This 
baryon asymmetry measurement is consistent with zero 
though it has large uncertainties, and central values that 
are large in magnitude. Note that the A[] asymmetry is 
expected to be dominated by hadronic effects [57]. For 
this reason it is not considered further in this work. The 
analogous effect at the LHC, cr{pp A°Ai) > a{pp —>■ 
Af,X), is seen at high rapidity by CMS [5S], and in the 
preliminary results of LHCb [29] . 

The predictions made in [13] for CDF are updated to 
include mixed EW-QCD corrections in an approximate 
way. These corrections are found to be small in magni¬ 
tude, and CDF’s measurements are found to be in good 
agreement with the SM predictions. The charge asym¬ 
metry at 7 TeV measured by LHCb is found to be in good 
agreement with Standard Model. It is also predicted that 
the charge asymmetry at 13 & 14 TeV will be smaller 
than at 7 & 8 TeV. The SM asymmetry is predicted to 
be very small for DO, which is consistent with what was 
measured. Several BSM models, including a 100 GeV 
axigluon model, are ruled out by this combination of 
SM predictions and measurements. On the other hand, 
it is shown that the Zbb couplings can be modified to 
explain the anomalous bottom-quark forward-backward 
asymmetry at LEPl m, while being consistent with the 
results of CDF and LHCb. 

The rest of the paper is organized as follows. In 
Sec. |n| the setup for the SM calculation of the forward- 
backward asymmetry is given. Predictions for the asym¬ 
metries measured by CDF, LHCb, and DO (and for those 
to be measured by LHCb) are given in Sec. |HI[ Fol¬ 
lowing that. Sec. jlVj discusses the approximation used 
for the mixed EW-QCD corrections to the asymmetry, 
the charm-quark charge asymmetry, and implications for 
BSM scenarios. Finally, conclusions are given in Sec. |V| 


II. STANDARD MODEL CALCULATION 

The Standard Model contribution to the top-quark 
forward-backward asymmetry has been investigated ex¬ 
tensively [TTl IT2l IT^ l3TH4T] .[65] Both the forward- 
backward and the charge asymmetries can be defined as 
follows. 


Afb,c 


cr (V > 0) - o- (V < 0) 
n (V > 0) -f o-(V < 0) ’ 


( 1 ) 


where a is the cross section in a given bin, and the 
observable Y is used to determine whether an event is 
“forwards” or “backwards.” The partonic level observ¬ 
able most closely related to what CDF measures is the 
so-called rest frame forward-backward asymmetry, corre¬ 
sponding to V = Vb — 2/5, while DO’s measurement is 
closely related to the lab frame asymmetry, Y = yb- 
Here is the rapidity of the (anti-)bottom quark, 

y = ln((i? + pz)/{E — pz))l2, and z is the proton di¬ 
rection at the Tevatron. The observable for the charge 
asymmetry measured by LHCb is V = 2/b ~ 2/f ■ 

The SM contributes to the forward-backward asymme¬ 
try at various orders in perturbation theory, which can 
be written schematically as 

Afb = (2) 

+ a^Ni + a^aNi + af V2 -!-••• 
a'^Do + + Off Di -I- olaDi -!-••• 

At the energy scales relevant for hadron colliders, Eq. ([^ 
can expanded in powers of coupling constants 


. JVi a JVo Jvi 

Afb =as — ^ — 1^ — -I- a— 
Dq aj Do Do 

N 2 NiDi' 


(3) 


Eq.(§ is the definition of the forward-backward asymme¬ 
try we will use in our calculations, as is commonly done. 
The approach is the same as that of m- We apply the 
analytic formulas in the literature for the 0{as) [H] and 
0[oA jo?^ [34] terms to the case of the bottom asym¬ 
metry. The Cuhre integration routine from the Cuba 
library |42j is used for the leading order and virtual-|-soft 
corrections, and MadGraph5_aMC@NLO [43] is used for 
the hard radiation. Unless otherwise stated, the set of 
parton distribution functions (PDF) used for all calcu¬ 
lations is the NNPDF2.3QED NLO grid [23J US] with 
as {Mz) = 0.119 and mb = 4.75 GeV. The other nu¬ 
merical values used in this analysis are: Mz = 91.1876 
GeV, Tz = 2.4592 GeV, a{Mz) = 1/127.940, and 
sin^ 6w = = 0.23126, which were taken from [46] . 

Cuts and binning are discussed in Sec. jlllj 

The tree-level EW contribution to the total cross sec¬ 
tion is included in the denominator of the computation of 
Afb.c when the invariant mass bin contains the Z-pole. 
Even though this term is formally higher-order than what 
we are considering, it is enhanced enough near « Mz 
to warrant inclusion. The numerically small flavor exci¬ 
tation contribution to the asymmetry, qg —/ bbq, is ne¬ 
glected in all cases except for the inclusive DO asymme¬ 
try. Similarly, t-channel VF-exchange is neglected for all 
bottom asymmetries, but is analyzed in the context of a 
charm-quark asymmetry. 

A formula for the N 2 term was not given in Ref. 
presumably due to its complicated nature, so we will 
drop the entire 0(as) term. As was done before N 2 
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was known, to compensate for this neglect of higher- 
order corrections we assign an uncertainty to the calcula¬ 
tion of 30% of the 0{as) contribution, originating from 
asDi Ki 0.3Z)o- Ref. [T7j finds that including Di, but 
not N 2 in the NLO QCD calculation of the top Aps de¬ 
creases the asymmetry by 25%, and that the full 0(a‘^) 
term increases the top asymmetry by 13%. 

In [T3] it was noted that the 0(a) contribution to 
is small compared to the 0{as) and 0[oA/o?^ terms, and 
their associated uncertainties. That statement is made 
quantitative in this work by computing the 0(a) contri¬ 
bution in an approximate way. The QED corrections 
are known to be in a one-to-one correspondence with 
the QCD asymmetry; one simply makes the following 
replacement for a given partonic channel: 

(4) 


However, the mass of the Z spoils this correspondence 
for the weak corrections. In this approximation, first 
one treats the Z as massless. Then, including real Z 
radiation, there is a one-to-one correspondence with the 
QCD asymmetry, 

12q; {T^ — 2QgS^) (T^ — 2QbS^) (5) 

This is what’s done in e.g. |35j . for the top asymmetry 
as 2mt ^ Mz- However, since [13] showed that the Z- 
pole is important for the bottom asymmetry, this result 
is then multiplied by a correction factor that attempts to 
account for the resonance structure of the Z^ 


{xf{x) f Re (1 - 


JdM^r {xf{x)) 


( 6 ) 


In Eq. 0, /r| = — iTzMz, and xf{x) is x times the 

PDE set for a given light quark flavor with x = 

The integration is over the range of for a given bin. 
The form of this correction factor is justified a posteri¬ 
ori due to the smallness of the O (a) terms; a different 
functional form won’t affect the result for the total asym¬ 


metry very much. This is discussed further in Sec. IV 


Note that an exact computation of the 0{a) contribu¬ 
tion to the ti asymmetry was made by [34]. Similarly, 
the 0{a) contribution to the bb asymmetry can be ex¬ 
tracted from m- In addition, shortly after this work 
was made public, a dedicated, exact computation for the 
0{a) contribution to the was given in Ref. [15]. 


III. RESULTS 

In this section, parton-level SM predictions are given 
for the forward-backward and charge asymmetry. The 
cuts and binning are tailored to match the experimental 
analyses. In the tables below, the superscript uncertainty 


comes from running the renormalization and factoriza¬ 
tion scales, fiR = from the chosen central value 

of Mz up to 2Mz. Similarly, the subscript uncertainty is 
due running /i down to Mz 12,. Sometimes there is a par¬ 
tial cancellation between the scale uncertainty from the 
tree level EW contribution and the NLO QCD contribu¬ 
tion. The QCD piece gets smaller as the scale increases 
because it is proportional to Og, whereas the EW piece 
gets bigger as /i increases because it scales as l/og. The 
additional, symmetric uncertainty in the last column is 
due to the neglect of higher-order QCD contributions. 
As previously stated, this uncertainty is taken to be 30% 
of the pure QCD asymmetry, as is typically done. 

As a check on the calculation, we have com¬ 
puted the NLO QCD contribution to the top- 
quark forward-backward asymmetry. We find 
A%p = 7.34%, Afg {Mtt< 450 GeV) = 5.36%, 

and A*pg [Ain > 450 GeV) = 10.59%, in good agreement 
with other determinations. 

In what follows, the rest-frame and lab-frame forward- 
backward asymmetries are denoted by A^pg and App 
respectively. 


A. CDF High Mass Analysis 

Results for the bottom-quark forward-backward asym¬ 
metry for the CDF high mass analysis [18j are given in 
Table [l] in percent. A cut has been placed on the ra¬ 
pidity of the bottom quarks, ly^jl < 1.1. No cuts have 
been made on the hard radiation. The results for two of 
three bins are in good agreement with CDF’s measure¬ 
ments. In the 225 — 325 GeV invariant mass bin, the SM 
prediction is just 1.3 ct above CDF’s result. 

The mixed EW-QCD corrections, 0{a), decrease the 
asymmetry. On the other hand, tree-level EW contribu¬ 
tion, which was neglected in the evaluation of the high- 
mass asymmetry in [18j . is positive though smaller in 
magnitude than the mixed EW-QCD corrections. The 
uncertainty due to the neglect of higher-order QCD terms 
is larger than the magnitude of the sum of the two EW 
contributions, and the scale uncertainty is comparable in 
magnitude (or larger than) the total EW contribution. 

The astute reader will notice that the predictions for 
the 0{as) asymmetry in Tab. |l] are slightly larger than 
those in Table I of [18], which used MSTW2008 NLO 
PDFs [15]. Indeed, though the two calculations agree 
within the scale uncertainty, the choice of PDF set seems 
to make a small difference in the prediction for the central 
value of the asymmetry. [66] For example, using Hollik 
and Pagani’s Eq. (7) [M], we find using NNPDF2.3QED 
NLO PDFs that the O contributions to A*^g 

are 2.5 • 10“^ and 3.6 • 10“^ in the uu and dd channels 
respectively. Using MSTW2008 NLO PDFs, we instead 
find 0{p?/cP'^uu = 1.8T0“^ and O^o? j = 3.2-10“"^, 
slightly smaller than the previous evaluation. Hollik and 
Pagani themselves find 2.3 • 10“^ and 3 • 10“^ in the uu 
and dd channels respectively [M] using MRST2001 LO 
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_Bin_ 0{a^lal) Ojas) Oja) A^b[%] 

150 < Mi,-,/GeV < 225 0.151^1 2.43;°;g® -0.15l°:°o 2.43 ± 0.73 ^^2 

225 < M,-,/GeV < 325 0.20in5 ^■72+^:li -O.SltU'.gl 4.61 ± 1.38 

325 < M,JGeV 0.281^6 8.99;g:^^ -0-57tg °^ 8.70 ± 2.61 


TABLE I: A“, 


in percent broken down by the various contributions and into different bins for the CDF high mass 
analysis. Here |?/^ < 1.1, and no cuts have been placed on the hard radiation. 


PDFs [SD]. 

B. CDF Low Mass Analysis 

Results for the bottom-quark forward-backward asym¬ 
metry for the CDF low mass analysis m are given in Ta- 
blelTT|in percent. The cuts PTb,b > 15 GeV and > 35 
GeV were placed on all of the bins. In Ref. m, only the 
Pt cut had been implemented. m The cut on the rapidity 
of the bottom quarks, ^ 1) is the same as in [T3]. In 
this work, two cuts were placed on hard radiation instead 
of one. In particular, the opening angle, A(j), between the 
bottom and anti-bottom in the plane transverse to the 
beam line is required to be greater than 2.8 radians. [55] 
In addition, to match CDF’s cuts, the bottoms are also 
required to satisfy, 

\PTb-PTb\ < 0.6max{pT&, Ptl}- (7) 

The cut in Eq. Q was not included in the analysis of 
Ref. HI]. 

No scale uncertainties are given for the 0{a) terms as 
they are smaller than 10“"^, which is due to the fact that 
the mixed EW-QCD corrections are themselves rather 
small. The 0{a‘^/a‘^) terms are all either in good agree¬ 
ment with the findings of m or are slightly higher. 
However, the bins that are slightly higher are exactly 
the bins that we expect to be affected by now having 
both a Pt and an cut. The affected bins are the 
35 < Mfjj/GeV < 75 bin, all three rapidity bins, and the 
35 < Mfjj/GeV bin. For the 0{as) terms, the cuts on 
the hard radiation increase the asymmetry by 8 — 27%, 
depending on the bin. As expected, without cuts on hard 
radiation, the asymmetries given here and in Ref. m are 
in good agreement. 


C. LHCb 


Results for the bottom-quark charge asymmetry for 
the LHCb 7 TeV analysis [33] are given in Table HI 


The cuts, 2 < j/f, f, < 4 and > 20 GeV, have 

been applied. In addition, the bottoms are required to 
have an opening angle in plane transverse to the beam 
line satisfying Acj) > 2.6 rad. The prediction in each 
bin for y/s = 7 TeV is in good agreement with mea¬ 
sured value. Without the inclusion of the tree-level EW 
contribution to the total cross section, the prediction 


for (75 < Alf^b/GeY < 105) increases from 3.50% to 
3.81%. 

Predictions have also been made for the bottom charge 
asymmetry at LHCb for ^/s = 8, 13, and 14 TeV. There 
are two differences between these calculations and the 7 
TeV analysis. First, the rapidity cut has been changed 
to 2.2 < Ub b ^ 4.2, and second, a low pt bin has 
been included. Results for the SM charge asy mme try 
at = 8, 13, and 14 TeV are given in Tables IV ^ 
and |VI| respectively. Due to the increase in gluon fusion 
initiated bottom production at higher s (smaller x for a 
given the asymmetry is predicted to be smaller at 

13 & 14 TeV than it is at 7 & 8 TeV. 

Unlike the case of the CDF low mass analysis, the Acj) 
cut makes very little difference for the LHCb analysis. It 
only increases the O (a^) asymmetry by 1—2%. The scale 
uncertainty is artificially small in the > 105 GeV bin 
due to a partial cancellation between the 0{a'^/a1) and 
0{as) terms. In fact, a peculiar feature of this cancella¬ 
tion is that for ^/s = 7 and 13 TeV, the asymmetry is the 
smallest in the Mf,i > 105 GeV bin for p = Mz rather 
than p = 2Mz or M^/2.[55] 


D. DO 

The same perturbative calculation used for the GDF 
and LHGb results is again used for the DO analysis. 
We are mindful that because DO measured an exclusive 
hadronic final state, a perturbative calculation may not 
necessarily be relevant. However, DO finds the that rms 
width of the distribution of {pb — Vb) is 0.11 |3D]; see [3T] 
for more information. This suggests that hadronization 
does not significantly affect Aps in meson produc¬ 
tion. 

The following cuts are made: 0.1 < \yb\ < 2.1 to match 
DO’s analysis, and \yi\ < 2.1 to simplify the calculation. 
Note that DO uses pseudorapidity in its analysis, whereas 
rapidity is used in this calculation. 

We find that the O (as) contribution to the inclusive 
asymmetry is 9.8±0.3 ^q ] • 10“®. The tree-level EW con¬ 
tribution to the asymmetry is 4.8^g g • 10“®. Since this 
asymmetry is so small, the flavor excitation piece is also 
considered. A quick calculation using MadGraph |i5] 
gives 0{as)gg = 1.1 ± 0.7(stat.) • 10“"^, the same order 
of magnitude as the qq initiated asymmetry. Gombin- 
ing these results yields our final prediction for the in¬ 
clusive asymmetry, App = (2.1 ± 0.8) • 10“"^. A very 
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Bin_ 0(o?lc?s) Ojas) Oja) A%[%] 


35 < Mfcs/GeV < 75 

u.uu_o.oo 

u.zu+o.oi 

-0.01 

0.19 ± 0.06 

75 < M^JGeV < 95 

2.oil°:r, 

o.52;°°32 

-0.05 

2.49 ± 0.16 

95 < MpJGeV < 130 

0.56lg:[I 

0.89;°°2" 

-0.01 

1.44 ± 0.27 

130 < M^i/GeV 

0.151^4 

9 1 1 -0-07 
^•-*--*-+0.08 

-0.13 

2.14 ± 0.63 

0.0 < \Ayp,\ < 0.5 

0.05lg°l 

f) 1 o-O.Ol 

U.lO+0.01 

-0.00 

0.18 ± 0.04 

0.5 < \Ayp,\ < 1.0 

nil +002 

U.ii_o,o2 

^•^y+0.02 

-0.01 

0.38 ± 0.09 

1.0 < \Ayp,\ < 2.0 

n 1O+0.03 

0 oq-0.02 

-0.02 

0.51 ±0.12 ini 

35 < M^-JGeV 

n nQ +0 02 

u.uy o 02 

^•^'->+0.02 

-0.01 

0.34 ±0.08 ini 


TABLE II: in percent broken down by the various contributions and into different bins for the CDF low mass 

analysis. No scale uncertainties are given for the 0{a) terms as they are smaller than 10“^. See the text for a 

discussion of the cuts used in this analysis. 


Bin 


0{as) 


40 < M^-JGeV < 75 
75 < MpJGeY < 105 
105 < M^i/GeV 

0.001°°° 

2.481°:®° 

0.25l°'°^ 

u.11:0+0.04 

-L.UZ+0.08 

1 Kr)—0.06 

-L.9O+0.09 

0.46 ± 0.14 ;° °l 

3.50 ± 0.31 Ig-jg 
1.79 ±0.46 mi 


TABLE III: in percent for LHCb with ^/s = 7 TeV. Here 2 < ^ < 4, Erp^^ t, > 20 GeV, and Aip > 2.6 rad. 


small asymmetry is exactly what’s expected from a sam¬ 
ple of bottom quarks without cuts on pT or Mf,j. There 
is nothing to suppress the symmetric gluon fusion pro¬ 
duction process, so the asymmetry is diluted away to 
almost nothing. DO’s result for the inclusive asymmetry 
is Afb = [-0.24 ± 0.41(stat.) ± 0.19(syst.)]% [20], 
which is in good agreement with this calculation. 

Predictions are also made for the |?7(i3)| and pt{B) dis¬ 
tributions of the asymmetry measured by DO, which are 
given in Fig. Aps is given as a function of \r]{B)\ in 
the left panel and aipriB) in the right panel. Data from 
DO and their corresponding predictions [20j . which were 
made using MC@NLO-|-Herwig are shown in black and 
purple respectively. The SM predictions from this work 
are in red. The absolute value of the rapidity and the pr 
of the bottom-quark are used for \r]{B)\ and pt[B) re¬ 
spectively. The flavor excitation process is not included 
in these distributions. Only the uncertainty due to the 
neglect of higher-order terms is included; no scale uncer¬ 
tainty is calculated for these distributions. The SM cal¬ 
culation is consistent with DO’s measurements in all bins 
except for the 7 < pr/GeV < 9 bin, where the asymme¬ 
try is measured to be larger than what is predicted. 

As previously noted, there is good agreement between 
the SM predictions in this work and the DO measure¬ 
ments. On the other hand, the DO observations and 
their predictions from MC@NLO-|-Herwig differ at the 3 ct 
level, with the MC prediction being larger than what was 
measured the majority of the time. Note also that the DO 
baryon analysis also finds that MC@NLO-|-Herwig pre¬ 
dicts an asymmetry that is larger than what was mea¬ 
sured ps] . 


IV. DISCUSSION 

In what follows a discussion is given of the estima¬ 
tion of the mixed EW-QCD corrections, the charm-quark 
charge asymmetry, and the implications of this work for 
BSM scenarios. 


A. Estimate of Mixed EW-QCD Corrections 

Inspecting Tables[l|and[n|it is seen that the 0{a) terms 
make small contributions to the bottom Aps relative to 
the 0{as) and 0[pA terms, and their associated un¬ 
certainties. This smallness justifies the functional form 
of the approximation used in Eq. § a posteriori, simply 
because changing the functional form won’t make much 
of a difference in the total prediction for the asymmetry. 
Case in point, consider the 95 < M^^/GeV < 130 bin 
in Table |ll| where there is a partial cancellation between 
the QED and the weak contribution, with (!I(q!)qed ~ 
—0.05% and ©(ajweak ~ -1-0.04%. Now consider what 
would happen if the following form for the weight was 
used, 

with s being the integration variable. [70| Such terms 
arise from interference between tree-level gluon exchange 
and box diagrams containing one Z and one gluon, see 
e.g. m- Replacing Mz with pz and taking the real part, 
the effect of the log is to change the sign of C>(a)weak, 
so that 0{a) « -0.05% - 0.04% = -0.09%.|7I] Con¬ 
sidering the 0{a) term by itself, this looks like a large 
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_Bin_ 0{o?lo?s) Ojas) A’S[%] 

10 < PTb. b/GeV < 20 O.OOtg'.^ 0.19 ± 0.06 

40 < M,5/GeV < 75 0.00t°;°° 0.48;n3 0.48 ± 0.14 
75 < M^i/GeV < 105 2.56^0.51 O.TS^Ht 3.31 ± 0.22 
105 < M,-JGeY 0.27^°;°^ 1.54;g:°g 1-81 ± 0-46 

TABLE IV: in percent for LHCb with ^/s = 8 TeV. Here 2.2 < yi, i < 4.2, Erp^ b > 20 GeV, and A(/) > 2.6 rad. 

_Bin_ 0{o?lal) Ojas) Ag[%] 

10<PTi,,5/GeV<20 O.OOtnS 0.10;^! 0.10 ± 0.03 
40 < M,b/GeV < 75 0.00t°:°° 0.177^2 0.17 ± 0.05 ;° 02 
75 < M,b/GeV < 105 1.47t° j® 0.44;HI 1-91 ± 0.13 
105 < M,-JGeY O.Mtg.g^ O.TO^Hs 0-84 ± 0.21 +°-°° 


TABLE V: Same as Table IV but for = 13 TeV. 


effect, almost an order of magnitude increase. However, 
the sum of all the contributions to the asymmetry in the 
95 < Mf,b/GeV < 130 bin only changes from 1.44% to 
1.34%. This change is smaller than the scale uncertainty, 
and it is also smaller than the uncertainty due to the 
neglect of higher-order QCD terms. 

B. Charm-Quark Charge Asymmetry 

There is also some interest in measuring the charge 
asymmetry in charm-quark production [SI IS] , and LHCb 
may well have the charm-tagging capabilities to do so. 


A full study of the SM contribution to the charm asym¬ 
metry is left for future work. Instead, in this work, we 
consider a previously neglected, tree-level contribution 
to the heavy quark charge asymmetry due to t-channel 
W exchange. For the top and bottom asymmetries, this 
contribution is rightfully neglected because it is highly 
Cabibbo-Kobayashi-Maskawa (CKM) (or bottom PDF) 
suppressed. However, there is less suppression for the 
charm-quark asymmetry because it is only a first to sec¬ 
ond generation transition. To the best of our knowl¬ 
edge, the asymmetric piece of the interference between 
s-channel gluon exchange and t-channel W exchange has 
not previously been given in the literature. 


_ asa \VQq\^ Cp tt/? c (3 -I- -|- (l - -I- 4 p^) -|- 8771“^ 

d cos 6 sin^ 9w Nc 8s (l + 2 p 2 _ 2m'^f - 


where = ttiqIs, = M'^/s, (3 = \/l — Am?, c = 
P cos 0, and V is the CKM matrix. The square of the 
t-channel diagram is suppressed by (a/as)(|14dP/s^) ~ 
1% relative to the interference term. The amplitude for 
ds (and sd) initiated cc production has less CKM sup¬ 
pression than the dd process. However, the ds + sd ini¬ 
tiated t-channel W exchange does not interfere with s- 
channel gluon exchange, so its effect on the asymmetry 
should also be small compared to the inference contribu¬ 
tion. 

Potential collinear singularities are regulated by the 
mass of the W, so based on counting powers of cou¬ 
pling constants, the size of the asymmetry due to W 
exchange should be suppressed by about an order of mag¬ 
nitude compared to the NLO QCD contribution to the 
asymmetry. We find the O (osO;) contribution to 
is -0.4% for 350 < Mcg/GeV < 950 and \yc,c\ < 1.84 


with Mw = 80.385 GeV and rric = GeV. Manohar 
and Trott find the NLO QCD contribution to is 
6.7% for 350 < M^g/CeV < 650 and 18% for 650 < 
Mcg/GeV < 950 [12]. This is consistent with the naive 
scaling estimate; t-channel W exchange is small, but not 
negligible for the charm asymmetry. On the other hand, 
we find the O {olsO) contribution to APpp for the same 
invariant mass range to be —6.10“®, which is completely 
negligible. 


C. Implications for BSM Scenarios 

Two different BSM scenarios are investigated in this 
work. First, as was done in [T3|, the effect on APpp of 
BSM physics models proposed for is considered. 

However, this analysis comes with the caveat that since 
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_Bin_ 0{o?lo?s) Ojas) A’S[%] 

10 < PTb,b/GeV < 20 0.00tg;|;g O.OS^Hi 0.08 ± 0.02 
40 < M,5/GeV < 75 0.00t°;°° 0.21;^^ 0.21 ± 0.06 
75 < M^i/GeV < 105 0.37;n3 1-72 ± 0.11 lol® 

105 < M,^/GeV 0.13tg:°^ 0.84;g:°6 0-97 ± 0.25 7°;°^ 


TABLE VI: Same as Table IV but for = 14 TeV. 


DO data 
SM (this work) 
MC@NLO+Herwig 


0.5 


1.0 

I?7(B)| 


1.5 



Pr(B) [GeV] 


FIG. 1: ApB vs. \r]{B)\ (left) and Pt{B) (right). In both cases, data from DO and their corresponding 
predictions [2D], which were made using MC@NLO+Herwig, are shown in black and purple respectively. The SM 

predictions from this work are in red. 


the anomalous top-quark forward-backward asymmetry 
has been resolved, none of the following models are nec¬ 
essarily viable anymore. Second, constraints on modified 
Zbb couplings are derived from measurements of c 
near the Z-pole by CDF and LHCb. A comparison of 
these bounds with the analogous results from LEPl |3D] 
is given as well. 

The forward-backward asymmetry as a function of in¬ 
variant mass is plotted in Fig. CDF’s low mass mea¬ 
surement m are in the top panel. The LHCb 7 TeV data 
is on the bottom left, and the CDF high mass data is on 
the bottom right. The SM predictions from this work 
are shown in red. Plotted in blue in Fig. is a predic¬ 
tion from the axigluon model |51j . The parameters used 
in the plot are, Mg> = 100 GeV, Pg' = Mc/lO, and 
Pa = 0.476. These parameters were taken from Ref. [52] . 
which throughly investigate the bounds on the axigluon 
models that were relevant as BSM explanations of 
The other parameter choices in Ref. [52] do not cause 
significant deviations from the SM, as can be seen in Fig. 
1 of (TD) . The measurements by CDF [TD| combined with 


the SM predictions in this work disfavors the 100 GeV 
axigluon. 

Furthermore, both of the benchmark points from |13j 
for the scalar weak doublet model |53| are also disfa¬ 
vored. This can be seen in Fig. where CDF’s mea¬ 
surements m are in black, and the SM predictions from 
Table |ll| are shown in red. In addition, the 100 and 150 
GeV flavor octet, electroweak triplet vectors of [M] IDS] 
are also disfavored, but 250 GeV vectors are consistent 
with this analysis. Predictions for the bottom asymme¬ 
try due to a 105 GeV scalar weak doublet and a flavor 
octet of 150 GeV EW triplet vectors are plotted in Fig.[^ 
in brown and green respectively. The BSM contributions 
to A^pb are computed using MadGraph [iS], and a sta¬ 
tistical uncertainty of 0.3% is included in addition to the 
SM uncertainties. 

The bb forward-backward asymmetry at LEPl, App\ 
was measured to be 2.3a below the SM prediction |30[ . 
This deviation can be explained by modifying the Zbb 
couplings as follows, 


C D 


e 

swcw 


{{Ti! 


Qb^W + SgbL) Pl + {-QbSw + ^9bR) Pr) b, 


( 10 ) 


where deviations from the SM are parameterized by 5gbL,R- A constraint on these modifications comes from 
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FIG. 2: Top: CDF’s low mass measurements of [TH]. Bottom left: LHCb’s 7 TeV measurements of [53] 
Bottom right: CDF’s high mass measurements of A^pg [T3|. The experimental data is shown in black. In all three 
plots, the SM predictions from this work are shown in red, and plotted in blue are predictions for the bottom 

asymmetry due to a 100 GeV axigluon. 


the ratio of the partial width Z —66 to the inclusive 
hadronic width of the Z at LEPl, which is consis¬ 
tent with the Standard Model prediction [30]. A two 
parameter fit of Sgi,L,R to Apg'^ and i?b is made using 
the theoretical and experimental values in [46] . Similar 
fits have been performed in the past, see e.g. [56]. The 
regions favored by the fit at the Itr and 2a levels are 
given in Fig. [^ in blue and orange respectively. CDF and 
LHCb have made measurements of A^g and near the 

Z-pole, which also constrains the parameters 6gbL,R- In 
Fig. 4a (left panel), the darker green and darker yellow 
regions correspond to values of SgbL,R that are consistent 
with both the CDF and the LHCb measurements at the 
Ict and 2a levels respectively. The lighter green region 


is allowed at Icr by CDF and 2a by LHCb. Lastly, the 
lighter yellow region is allowed by CDF at 2cr. Fig. |4b| 
shows a zoomed in version of |4a[ centered on the region 
allowed by LEPl. 

Fig. 1^ shows that the Zbb couplings can be modified 
to explain the anomalously low A^p'g while being con¬ 
sistent with the 66 asymmetry measurements at hadron 
colliders. This result is perhaps not so surprising since 
the measurements by CDF and LHCb have fairly large 
uncertainties. However, given the far larger amount of 
data expected during Run-2 of the LHC, it may possible 
for LHCb to constrain the parameter space for possible 
explanations of A^p'g . The CDF measurement favors 
slightly smaller values of SgbL for a given value of SgbR 
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CDF low mass data 
SM (this work) 

150 GeV Flavor Octet, EW Triplet 




50 


100 150 200 250 


300 


MbblGeV] 


MbblGeV] 


FIG. 3: CDF’s measurements of [19] are plotted in black, and the SM predictions from this work are shown in 
red. Plotted in brown and green respectively are predictions for the bottom asymmetry due to a 105 GeV scalar 
weak doublet and a flavor octet of 150 GeV EW triplet vectors. 



given in blue and orange respectively. Regions allowed by CDF and LHCb measurements of A’^pg and A’^ are shown 
in green and yellow. See the text for details on the parameter space allowed by CDF and LHCb. (right) Zoomed in 
version of the plot of the left, centered on the region allowed by LEPl. 


than the LHCb measurement does. However, the width V. CONCLUSIONS 

of the bands allowed by CDF and LHCb are about the 
same size. 

The preliminary results of CDF at both high and low 
invariant mass are consistent with the SM predictions 
made in this work and in Ref. m- The predictions 
of [13j were expanded on in this work to include the mixed 
EW-QCD corrections in an approximate way, which were 
found to be small in magnitude. 
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The charge asymmetry at 7 TeV measured by LHCb is 
found to be in good agreement with Standard Model. It 
is also predicted that the charge asymmetry at 13 & 14 
TeV will be smaller than at 7 & 8 TeV. In addition, it was 
shown that t-channel W exchange makes a non-negligible 
contribution to the charm-quark charge asymmetry. 

Both the preliminary results of CDF at low mass and 
LHCb results at 7 TeV include a measured asymmetry 
in the bin containing the Z-pole that is larger than the 
asymmetry in the adjacent invariant mass bins, as pre¬ 
dicted in this work and m- 
DO’s result for ApB is consistent with zero, and with 
the prediction of a very small asymmetry made in this 
work. On the other hand, the prediction for the inclusive 
asymmetry made by DO using MC@NLO-|-Herwig is 3.3 
standard deviations above what was observed. 

Several BSM scenarios proposed for A*pg , including an 
100 GeV axigluon, are ruled out by this combination of 
SM predictions and measurements. On the other hand. 


it was shown that the Zbb couplings can be modihed to 
explain the 2.3cr anomaly in Apg'^ at LEPl while be¬ 
ing consistent with the bb asymmetry measurements at 
hadron colliders. 
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